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Abstract 



Corrections to Einstein's equations that become important at small curvatures are 
considered. The field equations are derived using a Palatini variation in which the 
connection and metric are varied independently. In contrast to the Einstein-Hilbert 
variation, which yields fourth order equations, the Palatini approach produces second 
order equations in the metric. The Lagrangian L(R) = R — a 2 /R is examined and it is 
shown that it leads to equations whose solutions approach a de Sitter universe at late 
times. Thus, the inclusion of 1/R curvature terms in the gravitational action offers an 
alternative explanation for the cosmological acceleration. 
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Introduction 



One of the most interesting aspects of modern cosmology concerns the acceleration of 
the cosmological expansion. Recent supernovae [1, 2, 3, 4] and CMBR [5, 6, 7, 8, 9, 10] 
observations indicate that the expansion of the universe is accelerating, contrary to 
previous expectations. 

Most attempts to explain this acceleration involve the introduction of dark energy, 
as a source of the Einstein field equations. The nature of the dark energy is unknown 
but it behaves like a fluid with a large negative pressure. One possible candidate for the 
dark energy is a very small cosmological constant. Another, possibly related, problem 
involves the existence of dark matter. Observations of spiral galaxies, elliptical galaxies 
and galactic clusters indicates that these objects contain a large amount of dark matter. 
The difference between dark energy and dark matter is that dark matter clusters with 
the visible matter and dark energy is more or less uniformly spread throughout the 
universe. 

Of course, one possibility is that we do not understand gravity on these large scales. 
Since dark energy and dark matter are needed to explain phenomena in regions of low 
curvature we can attempt to modify Einstein's theory by adding corrections that become 
important when the curvature is small (see [13, 14, 15, 16] for other approaches that 
involve modifications of Einstein's theory). Recently two attempts [11, 12] to explain 
the cosmic acceleration along these lines have been made. They involve adding a term 
proportional to 1/R to the Einstein-Hilbert action and varying the action with respect 
to the metric (see [17, 18] for other papers that consider non-polynomial terms in the 
action). This approach leads to complicated fourth order equations that can be simplified 
by performing a canonical transformation and introducing a fictitious scalar field. It 
was shown in these papers that the modified field equations can produce the observed 
cosmological acceleration without the need of dark energy. 

In this paper I also consider a correction to the action that is proportional to 1/R, 
but I use the Palatini variational principle to derive the field equations. In this approach 
g^y and are taken as independent variables. In Einstein gravity the variation with 
respect to Y™ v gives the usual relationship between r° and the metric, i.e. Y^ v is the 
Christoffel symbol associated with the metric g^. The variation with respect to g^ gives 
RfiviX) — lg^ u R{T) = —kT^ v , where n = 8nG. Thus, the Palatini variation is equivalent 
to the Einstein-Hilbert variation in Einstein's theory. This is not the case however for 
other Lagrangians. In fact, it has been shown [19, 20] that the Palatini variation gives 
the usual vacuum Einstein equations for generic Lagrangians of the form L(R). This is 
to be contrasted to the purely metric variation that produces fourth order equations. If 
matter is included the Palatini variation still produces second order equations, but the 
are no longer identical to Einstein's equations. In this paper I show that the Palatini 
variation of the Lagrangian L(R) = R — a 2 /R, where a is a constant, leads to field 
equations that give an accelerating universe at late times. 
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The Field Equations 



The field equations follow from the variation of the action 

V9d 4 x (1) 



/ 
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where 



T>a a -pa ft r« _j_ fa -pa pA -pa /q\ 

Run = R^iiav R = Q^Rfiv, K = 87rG and L M is the matter Lagrangian. Here we consider 
a Palatini variation of the action, which treats g^ v and F" as independent variables. 
Varying the action with respect to g^ u gives 

L' (R)R^ - 1 -L{R)g iLV = -kT^ (3) 
where is the energy-momentum tensor and is given by 

Varying the action with respect to F^ v gives 

V a [L'y/gg^\ - ^V^L'Vgg^K - \v*[L'y/9Sr\% = 0- (5) 
By contracting over a and \i it is easy to see that this is equivalent to 

V a [L'(i?)^H=0. (6) 

This equation can be solved for the connection using a similar approach to that used in 
general relativity. Alternatively, one can define a metric = L </ M „ and it is easy to 
see that equation (6) implies that the connection is the Christoffel symbol with respect 
to the metric h^ v . A conformal transformation back to the metric g^ u gives (see [21] for 
the details on conformal transformations) 

= | ° v } + 2^[2*W - 9^9 a %L'\ (7) 

where the first term is the Christoffel symbol with respect to the metric g^ u . At first 
sight it might appear that this does not really define the connection since L contains 
derivatives of the connection. However, if we contract the field equation (3) we get 

RL'(R) -2L(R) = -kT. (8) 

If this equation can be solved for R = R(T), as we will assume here, then the terms 
in (7) involving L (R) can be expressed as derivatives of T. Since T contains only the 
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metric and not its derivatives the connection will involve only first derivatives of the 
metric and the field equations will then be second order in the metric g^ u . 
The Ricci tensor and Ricci scalar are given by 

= RM - ^(L'r 2 V^L'V u L' + (L'r^V.L' + ^(L')-V DL ' ( 9 ) 

and 

R = R{g) + SiL'y'nL' - ^(L')- 2 V ,L'WL' (10) 

where R^g) is the usual expression for in terms of g^ v and R = g^ u R^ v . 

If T = then the solutions to equation (8) will be constants and h^ u will be equal 
to g^y times a (positive) constant. This implies that R^ u = R^ u (g), and R = R(g). 
Thus, in a vacuum the field equations will reduce to the Einstein field equations with a 
cosmological constant for a generic L(R) (see [19, 20] for the vacuum case). The field 
equations (3) can be written in the Einstein form 



R^{g) - ^R{g)g^ 



—kS, 



(11) 



with a modified source Sn V given by 



S, u = (L'yXv-l- g(L')- 2 V M L'V,L' - (L'rV^V.L' - \g^Y^ + \ (r ~ ^ ) 9^ 

(12) 

Now consider the Lagrangian 



a 



(13) 



where a is a positive constant with the same dimensions as R and the factor of 3 is 
introduced to simplify future equations. The field equations for this Lagrangian are 



1 



a 



3R 2 

Contracting the indices gives 



R 



/IV 



R 



a 

3R 



R 2 - kTR - a 2 
and the solution to this algebraic equation is 

1 



R 



kT ± Vk 2 T 2 + Aa 2 



(14) 
(15) 

(16) 



For large \T\ we expect the above to reduce to R = kT, which follows from the Einstein 
field equations. Thus, if T > we need to select the positive sign and if T < we need 
to select the negative sign. In a universe filled with an ideal fluid T = — (p — 3P), so 
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that T < if the dominant energy condition holds (i.e. p > and p > 3\P\ 
have 

R = \[kT - Vk 2 T 2 + 4a 2 } . 



Thus we 



(17) 



The vacuum solution is 



R = -a (18) 
0, the universe will approach a de Sitter spacetime and 



so that at late times, as T^ v 
the expansion of the universe will accelerate. 

From equations (14) and (16) we see that the field equations reduce to the Einstein 
equations if \kT\ » a. Thus, in a dust filled universe the evolution will be governed by 
the Einstein field equations at early times. Eventually the corrections to the equations 
of motion will become important and the universe will make a transition to a de Sitter 
universe at late times. To match the observations of the cosmological acceleration we 
must take a ~ 1CT 67 (eV) 2 ~ 1CT 53 m" 2 . Note that in a vacuum we get the Einstein 
field equations plus a small cosmological constant, so that this theory will pass all the 
solar system tests that general relativity has passed. It is also interesting to note that 
in a radiation dominated universe T = so that the dynamics is not governed by the 
Einstein's equations even at large curvature. As we will see below (see equation (24)) 
the equations of motion are the Einstein field equations with a cosmological constant 
and a modified Newton's constant. 

Now consider the evolution of a universe with metric 



ds 2 = -dt 2 + a(t) 2 [dx 2 + dy 2 + dz 2 } 
at late times, when a » kT. In this regime 



Rp, ^ R»M + jj^^uL' + ^jg^DL', 

R ~ -kT - a, 
2 



L ~ — a 
3 



1 - 



kT 



a 



L ~ 



1 + 



kT 
4a 



and the field equations are 



RfJLV — ^ OtQllV K 



3 T 5 

4 " ~ 16 ^ 



(19) 

(20) 

(21) 
(22) 
(23) 

(24) 



The matter in the present universe can be approximated by dust with T = po/o, 3 , where 
Po is a constant. The nonvanishing components of the Ricci tensor are 



a 8aa 3 



a / a\ ' 
--3 - 

a \a/ 



(25) 
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and 



Rij 



ad + 2d + 



. 2 3/tpo (o- & 2 



8aa Xa^a 2 , 



(26) 



At late times the universe will almost be in a de Sitter phase and we can take 

a(t) = e Ht + b(t) (27) 
where \b(t)\ « e Ht and a = 12H 2 . To lowest order in b the field equations are 

b-H 2 b=-^e- 2Ht (28) 
12 v ; 

and 



b + AHb-5H 2 b=^e- 2Ht . 

4 

Subtracting these two equations gives the first order equation 

b-Hb = ^e- 2Ht , 
12H 

and the particular solution to this equation is 

b(t) = 



K P0 -2Ht 



3QH 2 



(29) 



(30) 



(31) 



Thus, at late times the universe approaches a de Sitter spacetime exponentially fast. 
This behaviour is analogous to the cosmic no hair theorem for fourth order gravity 
discussed by Kluske and Schmidt [22]. 



Conclusion 

Using a Palatini variation the field equations for a nonlinear gravitational Lagrangian 
coupled to matter were found. The vacuum field equations are the Einstein equations 
with a cosmological constant. Thus, at late times as — > our universe will approach 
a de Sitter spacetime. The inclusion of matter gives field equations that differ Einstein's 
equations. Using these equations it was shown that the approach to de Sitter space is 
exponentially fast when the 1/R term dominates. Thus, the inclusion of non-polynomial 
curvature terms in the gravitational action offers an alternative explanation for the 
cosmological acceleration. 
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